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Abstract 

An equation of state (EOS) for uniform asymmetric nuclear matter (ANM) is 
constructed at zero and finite temperatures by the variational method start- 
ing from the nuclear Hamiltonian that is composed of the Argonne vl8 and 
Urbana IX potentials. At zero temperature, the two-body energy is calcu- 
lated with the Jastrow wave function in the two-body cluster approximation 
which is supplemented by Mayer's condition and the healing-distance condi- 
tion so as to reproduce the result by Akmal, Pandharipande and Ravenhall. 
The energy caused by the three-body force is treated somewhat phenomeno- 
logically so that the total energy reproduces the empirical saturation con- 
ditions. The masses and radii of neutron stars obtained with the EOS are 
consistent with recent observational data. At finite temperatures, thermo- 
dynamic quantities such as free energy, internal energy, entropy, pressure 
and chemical potentials are calculated with an extension of the method by 
Schmidt and Pandharipande. The validity of the frozen-correlation approx- 
imation employed in this work is confirmed as compared with the result of 
the fully minimized calculation. The quadratic proton-fraction-dependence of 
the energy of ANM is confirmed at zero temperature, whereas the free energy 
of ANM deviates from the quadratic proton-fraction-dependence markedly 
at finite temperatures. The obtained EOS of ANM will be an important 
ingredient of a new nuclear EOS for supernova numerical simulations. 
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1. Introduction 

The nuclear equation of state (EOS) plays important roles in the studies 
of high-energy astrophysical phenomena: The nuclear EOS at zero temper- 
ature governs the structure of cold neutron stars (NSs), whereas the EOS 
for finite temperatures is necessary for studies of core collapse supernovae 
(SNe), black hole formations, cooling of NSs and so on. At zero tempera- 
ture, a variety of nuclear EOSs, including those based on the realistic nuclear 
Hamiltonian, have been applied to NSs, and recent observational data on NSs 
impose severe constraints on the EOS [2|. On the other hand, the nuclear 
EOSs available for the numerical simulations of SN are limited, because it 
is a hard task to construct an EOS for the SN simulations: The SN-EOS 
table must provide thermodynamic quantities in an extremely wide range 
of densities, temperatures and proton fractions. Furthermore, the EOS of 
non-uniform nuclear matter as well as that of uniform matter is important 
in SN simulations. 

For the above reasons, until several years ago, there were only two kinds of 
nuclear EOSs available for SN simulations. One was constructed by Lattimer 
and Swesty jsj with use of a Skyrme-type interaction for uniform matter and 
a compressible liquid-drop model for non-uniform matter. The other was 
constructed by Shen et al. 01 with use of a relativistic mean field (RMF) 
theory for uniform matter and a Thomas-Fermi (TF) approximation for non- 
uniform matter. Recently, new SN-EOSs based on the RMF theory have been 
proposed: For the high-density region, the Shen-EOS was extended so as to 
take into account hyperon mixing ^] and the quark-hadron phase transition 
jcj. Furthermore, the RMF calculations with parameters different from those 
of the Shen-EOS have been applied to SN-EOSs 0, [sj. For the low-density 
region, the nuclear statistical equilibrium is taken into account 0, Isl , or the 
virial expansion is adopted (8|. 

These recent advances in the SN-EOSs are astonishing. However, in these 
EOSs, uniform matter is treated with phenomenological nuclear models: To 
our knowledge, there are no EOSs based on the nuclear Hamiltonian com- 
posed of bare nuclear forces. This situation motivates us to construct a new 
nuclear EOS for SN simulations based on a realistic nuclear Hamiltonian with 
use of the variational many-body theory. 
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In the standard variational method, the Jastrow wave function is as- 
sumed, and the expectation value of the Hamiltonian is cluster-expanded. 
Then, parts of the higher-order cluster terms are resummed by the Fermi Hy- 
pernetted Chain (FHNC) method T^. Akmal, Pandharipande and Raven- 



hall (APR) performed the FHNC calculation with the Argonne vl8 
(AVI 8) two-body potential 12| and the Urbana IX (UIX) three-body po- 
tential [13[; the obtained EOS has been referred to as one of the standard 
nuclear EOSs at zero temperature. In Ref. [l4|, the EOS of APR is extended 
to finite temperatures with the method proposed by Schmidt and Pandhari- 



pande (SP) 15|, |l6|. Furthermore, the quantum Monte Carlo (QMC) tech 



niques have been developed for nuclear matter at zero temperature Il7|,ll8| 



and, in Ref. [19|, free energy of hot nuclear matter is calculated with the 
method of SP in connection with the EOS by one of the QMC methods. How- 
ever, these variational methods are applied only to limited cases: symmetric 
nuclear matter (SNM) with the proton fraction x being x = 1/2 and pure 
neutron matter (PNM) with a; = 0. This is because the energy of asymmet- 
ric nuclear matter (ANM) with 0< x <l/2 is difficult to calculate by those 
variational methods. In fact, there are only few examples of variational cal- 
culations for ANM [i^, 21, 22 1. Furthermore, sophisticated numerical many- 



body calculations are, practically, not suitable for constructing an SN-EOS, 
because, including the non-uniform phase, hundreds of thousands of data on 
thermodynamic quantities must be completed in an SN-EOS table. 

In our project, we adopt a relatively simple variational method that can 
treat ANM, and have so far constructed a realistic nuclear EOS for SNM and 
PNM in Ref. 23|], hereafter referred to as paper I. In that paper, starting from 



the nuclear Hamiltonian with the AVIS and UIX potentials, we calculated the 
two-body energy in the two-body-cluster approximation, supplemented by 
appropriate constraints, and incorporated the three-body energy somewhat 
phenomenologically so that the total energy could reproduce the empirical 
saturation data. At finite temperatures, we employed the method of SP to 
obtain the free energy. 

In order to complete an SN-EOS, it is necessary to treat non-uniform 
matter as well as uniform matter. Since we focus on improving the treatment 
of uniform matter using the variational method, we simply adopt the TF 
calculation for non- uniform matter, following the method by Shen et al 
Here, we take into account the isolated atomic nuclei: Since a large amount 
of experimental data on the atomic nuclei are available, we fine-tuned the 
EOS so that the TF calculation for the atomic nuclei could reproduce the 
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gross features of their masses and radii in Ref. [2J], hereafter referred to 
as paper II. In that paper, we further constructed a non-uniform /9-stable 
EOS at zero temperature, which corresponds to the EOS of the NS crust, 
by the TF calculation. The obtained EOS is consistent with typical EOSs 



calculated for NS crusts [25|, |26 



Toward an SN-EOS, the next important task is to construct the EOS of 
uniform ANM with < a; < 1/2, which is the main object of this paper. 
Namely, we extend the variational method in papers I and II to obtain the 
EOS of uniform ANM at zero and finite temperatures. In Section 2, we per- 
form this extension for nuclear matter at zero temperature, and calculate the 
EOS of ANM. Then, we apply the EOS to NSs, and compare the NS solutions 
with recent observational data. In addition, we discuss the x-dependence of 
the energy of ANM quantitatively. In Section 3, we treat ANM at finite 
temperatures to calculate various thermodynamic quantities of ANM. We 
also examine the validity of the so-called frozen-correlation approximation 
adopted in our study: For this purpose, we perform the full minimization 
of the free energy in the present variational method. Furthermore, we dis- 
cuss the x-dependence of the free energy of ANM at finite temperatures. A 
summary is given in Section 4. 

2. Asymmetric nuclear matter at zero temperature 

2.1. Extension of the variational method to asymmetric nuclear matter 

In this section, we calculate the energy per nucleon E/N oi uniform ANM 
as a function of the total nucleon number density p and the proton fraction x. 
As in papers I and II, we start from the non-relativistic nuclear Hamiltonian 
H for nucleons interacting through realistic two-body and three-body 
forces, and decompose H into the two-body Hamiltonian H2 and the three- 
body Hamiltonian ifs as follows: 

H = H2 + H,. (1) 

The two-body Hamiltonian H2 is written as 

N k2 N 

^2 = -E^v? + E% (2) 

where m is the nucleon mass: For simplicity, we choose m in Eq. ([2]) as 
the mass of a neutron, as in paper II. For the two-body nuclear potential 
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Vij acting on an nucleon pair, we use the isoscalar part of the AV18 
potential as in papers I and II, namely, 

1 1 

E E [ycts{rij) + sVTt{rij)STij + sVsot{rij){Lij ■ s) 

t=Os=0 

+VqLts{rij) \Lijf + sVqsot{rij){Lij ■ s)^]Ptsij. (3) 

Here, S^ij is the tensor operator, Lij is the relative orbital angular momen- 
tum operator, and Ptsij is the spin-isospin projection operator projecting a 
two-nucleon state onto the spin-isospin eigenstates with t and s being the 
two-nucleon total isospin and spin, respectively. 
The three-body Hamiltonian H-^ is given by 

N 

H3= E V,,k. (4) 

i<j<k 

As the three-body nuclear potential Vijk-, we employ the UIX potential, which 
consists of the repulsive component V^^^, and the two-pion exchange compo- 
nent V^^^: 

Vr.k = V^, + VSl. (5) 

We first calculate the two-body energy E2/N without the three-body 
force. As in papers I and II, we use the two-body cluster approximation for 
the expectation value of H2 with the following Jastrow wave function: 

^F. (6) 

Here, is the degenerate Fermi-gas wave function at zero temperature 
and Sym is the symmetrizer with respect to the order of the factors in the 
products. The two-body correlation function fij is taken as 

h = EEE [fLir^J) + s/^,(r.,)5T., + sf^o,{r,,){L,, ■ s)] P,%. (7) 

t=0 /i s=0 

The Jastrow wave function of ANM in Eq. ([6]) differs from those of symmet- 
ric nuclear matter (SNM) and of pure neutron matter (PNM) in the following 
ways: i) The Fermi wave numbers of protons fcpp = (Svrpa;)^/^ and neutrons 
^Fn = [37rp(l — x)]^/^ included in are in general different from each other. 



^ = Sym 
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depending on x. ii) The correlation function fij depends on the third compo- 
nent of the isospin ^3: Even though the two-body interaction Vij is isoscalar, 
we expect that the two-body correlations for isospin-triplet three states (p-p, 
p-n, n-n) may differ from each other in, e.g., extremely neutron-rich ANM, 
due to an effect of the medium. In Eq. ([7]), the ts-dependence of fij is ex- 
pressed with the superscript /i, i.e., fi = (+, 0, — ) for ^3 = (1, 0, —1). Namely, 



the central, tensor and spin-orbit correlation functions 



), f^tivij) and 



depend not only on (t, s) but also on /i. In addition, the projection 
operator Pl'^-j projects the nucleon pair states onto the fi = (-|-,0,— ) 
states. Therefore, in the case of ANM, these 16 correlation functions are 
regarded as variational functions. 

With use of this trial wave function ip", we express the expectation value 
of H2 in the two-body cluster approximation to obtain E2/N: 



N 



El 
N 



t=0 fi s=0 Jo 



sF^,(r)VTt(r) 



r^dr[Fi:{r)Vctsir) 



m t=o fj. s=oJo 



qhts' 



r dr 



Cts\ 



n 2 



+48s 

On the right-hand side of Eq. 
energy given by 

Ei{p,x) ^ 



dr 



+ 8s 



df. 



dr 



dr 



^ qFts 



(8) 



(j8]), the first term Ei/N is the one-body kinetic 



X 



2m 



x] 



' '■ ^Fn 

2m 



(9) 



The second term represents the potential energy, where F^g(r), FJ^^{r), F^Q^^r), 



F^usif) and Fq^soi(^) are defined as 



Etsir) 



FsUr) 



■4 



Ftl 



-2A[fUr)fFl},,{r), 



qFtl'' 



(10) 

(11) 



(12) 
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F!:us{r) = 48.[/^,(r)]^Fp^,,(r) + {[/^,,(r)]V8.[/^,(r)]^}F, 



qFts('") 



+ ls[f^otir)?Fi:^Jr), (13) 



fUr)f^otir) + f If^M" + ^/T".(r)/s%(r)}F^'^p,i(r) 

+|[/s"o*W]'i^b%W- (14) 

Here, Fp^g(r), F^p^^{r) and -FbFts(^) various distribution functions in the 
case of the degenerate Fermi gas, the definitions and exphcit expressions of 
which are given in Appendix A. We note that these expressions are similar 
to those for SNM and PNM shown in paper I, though the x-dependence is 
introduced into the relevant functions with the superscript /i. 

Then, we minimize E2/N with respect to /ctsl'^); fTti''^) /soi(^) 
solving the Euler-Lagrange equations. Here, as in papers I and II, we impose 
two constraints in order to compensate effectively for the neglect of higher- 
order cluster terms. 



The first constraint is the extended Mayer's condition 23|, |27|, which is a 
kind of normalization condition, expressed as 

POO 

P mr)-F^Jr)]dr = 0. (15) 
Jo 

Here, F^^ir) is a two-body cluster approximation of the (t, s, /i)-projected 
radial distribution function shown in Eq. ( ITOj) . When Eq. ( !T5|) is summed 
over /i, this condition reduces to that adopted for SNM and PNM in papers 
I and II. 

The second constraint is the healing-distance condition expressed as 

fLir) = l. fUr)=0, /s%(r) = (r > n,), (16) 

which means that the correlation between two nucleons vanishes at the dis- 
tances r that are larger than the healing distance r^. Here, we assume that 
Th is common to the central, tensor and spin-orbit correlation functions, and 
proportional to the mean distance between nucleons, i.e., rh = a^ro with 
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ro = (3/47rp)-^/^. The coefficient is chosen to be = 1.76 so that E2/N 
reproduces the result of the FHNC calculation by APR which includes 
the higher-order cluster terms. As shown in paper I, the obtained E2/N 
with Oh = 1.76 for both SNM and PNM are in good agreement with those 
by APR, and, in this paper, we assume that Oh is independent of x. 

Next, we calculate the three-body energy E-^/N, based on the three-body 
Hamiltonian H^. Since there is some uncertainty in the three-body force, we 
express E-^/N somewhat phenomenologically as 

= ^( E + ^VSl]h + ip'e'^" [1 - (1 - 2x)^] . (17) 

i<j<k 

Here, the brackets with the subscript F represent the expectation value with 
the degenerate Fermi-gas wave function. The parameters a and /3 represent 
the effects of the dynamical correlations among nucleons and the relativistic 
correction ll|]. The last term on the right-hand side of Eq. ( !T7|) is a phe- 
nomenological correction term. This E^/N for ANM is a straightforward 
extension of those for SNM and PNM given in papers I and II, as explained 
below. 

The ffist term on the right-hand side of Eq. (IT7|) depends on x, due to 
the x-dependence of the degenerate Fermi-gas wave function for ANM. Since 
we have chosen in papers I and II to treat a and /3 as independent of x 
to obtain reasonable values oi E/N for both SNM and PNM, we make the 
same assumption for these parameters in this paper. In contrast, in paper 
I, we introduced a correction term corresponding to the second term on the 
right-hand side of Eq. (fT7|) in the case of SNM only, and set 7 = for PNM, 
following the treatment by APR. In this paper, therefore, we simply assume 
that the correction term goes to zero quadratically as x approaches zero, as 
explicitly shown in the last factor in Eq. fll7p : The x-dependence of the 
energy of ANM will be discussed in more detail later. 

Finally, the total energy per nucleon E/N is expressed as 

E Eo E'i , , 

The parameters a, /3, 7 and 5 in Eq. (JTTj) were determined in paper I so 
that the obtained E jN could reproduce the empirical saturation density 



-'^The validity of this treatment of the three-body force is examined for SNM and PNM 
inRef. [H. 
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Figure 1: Energies per nucleon, E/N, for various values of the proton fraction x as func- 
tions of the number density p. The energies obtained by APR for SNM and PNM are also 
shown. 

saturation energy Eq/N, incompressibility K and symmetry energy -Egym = 
E{pq,x = 0)/N — E{pq,x = 1/2) /N. Furthermore, in paper II, the values 
of these parameters were fine-tuned so that the Thomas-Fermi calculations 
for isolated atomic nuclei with the obtained E/N could reproduce the gross 
features of the experimental data in their masses and radii: For our purpose 
of constructing a reliable SN-EOS, this fine-tuning is an important process. 
Exphcitly, a = 0.43, /3 = -0.34, 7 = -1804MeVfm*^ and 6 = 14.62fm3 to 
obtain po = O.iefm'^, Eq/N = -16.09 MeV, K = 245 MeV and E^y^ = 30.0 
MeV.I 

The obtained E/N for ANM are shown in Fig. [T] Also shown are the 
results of APR. In the case of SNM (x = 0.5) at p < 0.32fm-^ the present 
E/N is in good agreement with that obtained by APR. At p > 0.32fm~^, 
however, our E/N becomes higher than that obtained by APR, due to the 7r° 



^We note that, due to the improvements in the numerical calculations, the values of the 
physical quantities shown in this paper are slightly different from those previously given 
in papers I and II . 
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condensation occurring in the study by APR. As x decreases, E /N increases 
monotonically. Then, in the case of PNM at x = 0, the present E/N is lower 
than that obtained by APR except for in the low-density region, because the 
symmetry energy obtained by APR is about 34 MeV, which is larger than 
the present result. 

Next, we apply the obtained E{p,x)/N to neutron stars (NSs). As in 
papers I and II, we treat NS matter as a charge-neutral, /3-stable mixture 
of nucleons, electrons and muons at zero temperature. For the NS crust, we 
employ the crust EOS calculated by the Thomas-Fermi method in paper II. 
We note that the obtained EOS of NS matter violates the causality, i.e., the 
sound velocity exceeds the speed of light, at densities higher than the critical 
density Pc = 0.88 fm~^, which is close to the corresponding value of the EOS 



by APR pjj. In this paper, therefore, we construct a modified causal EOS 
by following the method of Ref. ^29!], i.e., at densities higher than pc, we 
replace the pressure by P{e) = P{ec) + e — Cc, where e is the energy density 
and Ec is the energy density at pc- 

With the original and modified EOSs, we solve the TOV equation to 
obtain the masses and radii of NSs. Figure [2] shows the obtained mass-radius 
relation of NSs. It can be seen that the stable NS solutions with the original 
and modified EOSs are hardly distinguishable from each other, and both the 
maximum masses are 2.22Mq. Also shown in this figure are the mass of 
the heavy NS J1614-2230 (M = 1.97 ± 0.04Mq) Q and the observationally 
suggested region analyzed in Ref. j2|. The masses and radii of NSs with the 
present EOS are consistent with those of the observational data. 

2.2. Proton- fraction- dependence of the symmetry energy 

In this subsection, we discuss the proton- fraction-dependence of E/N and 
the symmetry energy -Egym in more detail. For this purpose, we define the 
following energy difference: 

EM^^-^^^^. (19) 

From this definition, it is obvious that Esym = Ei{po, 0). In many cases, this 
symmetry energy -Egym is regarded as the same as 5*1(^0) given by the original 
definition, i.e.. 



8dx^ N 



(20) 

x=l/2 
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Figure 2: Mass-radius relation of NSs with the present EOS. The filled circle represents 
the NS for which the central density is equal to the critical density pc: For smaller radii, 
the solid and dotted lines show the masses of NSs with the original EOS and the modified 
causal one, respectively. The horizontal band shows the mass of PSR J1614-2230 (1.97 ± 
O.O4M0) The shaded region is given in Ref. [5] based on the observational data: The 
dark and light gray regions correspond to la and 2a errors, respectively. 



These two quantities, -Esym and 5*1 (po), are identical if, as is usually assumed, 
E{p,x)/N increases quadratically with x from x = 1/2 (SNM) to x = 
(PNM), or equivalently, if Ei{p,x) is proportional to C = (1 — 2x)^. 

In order to examine this x-dependence of Ei{p,x) in the present calcula- 
tion, we show El at various densities as functions of ( in Fig. [31 It can be 
seen that, in the wide range of densities, Ei is nearly proportional to (. 

For a more quantitative examination, we try to evaluate E{p,x)/N by 
interpolating the energy per nucleon for PNM, E{p, 0)/N, and that for SNM, 
E{p, 1/2) /N, as 



£;int(p,x) E{p,l/2) 



N 



N 



+ 



E{p,Q) Eip,l/2) 



N 



N 



{l-2xf. (21) 



Then we show in Fig. Hlthe relative deviation of -Eint(p, x)/N from E{p, x)/N, 
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defined by 

/.p. s E,,,,{p,x)/N-E{p,x)/N 

= M^) ■ ^22) 

It can be seen that AE{p, x) is concave with respect to (, and < AE < 0.01 
at the densities shown in this figure. The main cause of this concaveness is 
the one-body kinetic energy Ei/N which is a convex function of The 
dotted curves in Fig. |4]will be explained in Appendix B. 

Figure [5] shows the density dependence of Si{p) and Ei{p,x). As can be 
seen in Fig. [5], Si{p) is shghtly lower than Ei(p, 0), because, in Fig. [3], Ei{p, x) 
is slightly convex with respect to (, and Si{p) corresponds to the gradient of 
Ei{p,x) at C = 0. The shaded region in Fig. [5] shows the experimental data 
on 5*1 (p) obtained from heavy-ion collisions |30|; the symmetry energy of the 
present EOS is consistent with those experimental data. 
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Figure 4: Relative deviations AE of the interpolated energies from the original ones as 
functions of ( (Type 1). The dotted curves show AE with the interpolation by Lagaris 
and Pandharipande, which is discussed in Appendix B (Type 2). 



3. Asymmetric nuclear matter at finite temperatures 

3.1. Extension of the variational method by Schmidt and Pandharipande 

In this subsection, we calculate the free energy per nucleon F/N of ANM 
as a function of the number density p, proton fraction x and temperature T. 
To obtain F/N for SNM and PNM, in paper I, we adopted the variational 
method by Schmidt and Pandharipande (SP) [l5|]. In this method, F/N is 
expressed with the effective mass of a nucleon m*, and then F/N is minimized 
with respect to m*, as explained below in more detail. The validity of this 



method was examined in Ref. 16|], and has been employed in some many 



body calculations of hot SNM and PNM [IJ, [191, |31[. Since our purpose is 
to treat ANM, we extend the method by SP to ANM, i.e., we distinguish 
the effective masses of a proton and a neutron m* {i =p, n), and treat them 
independently. 

We start from the following expression: 

l = ^-T^, (23) 
N N N' ^ ^ 

where Eq^/N is the approximate internal energy per nucleon and Sq/N is 
the approximate entropy per nucleon at temperature T. Here, the adjective 
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Figure 5: The density-dependent symmetry energy Si{p) defined by Eq. (PU)) as a function 
of the density p. Ei{p, x) are also shown for various values of x. The shaded region shows 
the experimental data on Si{p) (30| . 



" approximate" is attached in order to distinguish these quantities from those 
calculated with F/N through the thermodynamic relations. 

As in the case of zero temperature, we express Eq^/N as a sum of the 
two-body internal energy E2t/N and the three-body internal energy E^t/N: 



E, 



OT 



E2T -Est 



A^ 



A^ 



(24) 



At zero temperature, E2T/N reduces to E2/N given in Eq. (|8]), where the 
occupation probabilities of single-particle states in the Jastrow wave function 
Eq. ([6]) are expressed as nQi{k) = 6{k-pi — k) {i = p, n). At finite temperatures, 
we follow the method of SP, and construct E2T/N by replacing noi{k) in E2/N 
at zero temperature by the following averaged occupation probabilities ni{k): 



tiiik) 



1 + exp 



-1 



(25) 
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where ei{k) are the quasi-nucleon energies which are assumed to be as follows: 



e^{k) = (26) 

Here, m* [i = p, n) are the effective masses for protons and neutrons, re- 
spectively. In Eq. (l25l) . the values of /ioi are determined by the following 
normalization conditions: 

1 f°° 

= — / n,{k)k^dk, (27) 
Jo 

where pp and pn are the proton and neutron number densities, respectively. 
The explicit expression of E2t/N obtained in this method is similar to E2/N 
at zero temperature given by Eq. ([H]), but Ei/N, Fpj^(r), -FqFts(^) 
F^^Jr) are replaced by E,^/N, F^Jr;T), F^^Jr;T) and F,V(r;T), the 
explicit expressions of which are given in Appendix A. Following the method 
by SP, the two-body correlation functions fcts{r), f^^ir) and f^otij") are as- 
sumed to be the same as at zero temperature: We refer to this approximation 
as the frozen-correlation approximation (FCA), the validity of which will be 
discussed later. 

The three-body internal energy E^t/N in Eq. flM|) is assumed to be the 
same as E^/N at zero temperature in Eq. (fljl) . as in paper I, where the 
validity of this assumption is also discussed (23|. 

In the method of SP, the approximate entropy Sq/N in Eq. f l25]) is ex- 
pressed as in the case of a gas of quasi-nucleons whose masses are the effective 
mass m*. In this paper, we extend the method so as to treat ANM, i.e., 

§ = - E^ ^{[l-n,ik)]^n[l-n,ik)]+n^{k)\nn,{k)}edk. (28) 

i=p,n ^ P 

Then, F/N is expressed as an explicit function of m* and m*, and is mini- 
mized with respect to them. 

3.2. Free energy and related thermodynamic quantities 

The obtained free energies per nucleon of ANM are shown in Fig. E) 
Also shown are those obtained with the FHNC calculations for SNM and 



PNM given in Ref. [ij], which are extensions of the results of APR. As 



in the case at zero temperature, the present result is in good agreement 
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Figure 6: Free energies per nucleon at T = 20 MeV (a) and T = 30 MeV (b) for various 
values of x as functions of p. The free energies calculated with the FHNC method for 
SNM and PNM given in Ref. 



14| are also shown. 
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Figure 7: Internal energies per nucleon, Ej^/N, at T =10 and 30 MeV for various values 
of X as functions of p. The approximate internal energies per nucleon, Eq'y/N, are also 
shown. 

with that obtained with the FHNC method for SNM. As x decreases, F/N 
increases monotonically, and, in the case of PNM, our result is somewhat 
lower than that calculated with the FHNC method because of the larger 
symmetry energy in the case of APR, as pointed out in the last section at 
zero temperature. 

Figure [7] shows the internal energy per nucleon E-^/N of ANM derived 
from F/N through the thermodynamic relation. E^/N increases as x de- 
creases monotonically from SNM to PNM. Also shown is the approximate 
internal energy Eq^/N. It can be seen that E^/N is in good agreement with 
Eqt/N for various values of x, which implies that the present variational 
method is self-consistent at finite temperatures. In paper I, we confirmed 
this self-consistency for SNM and PNM: In this paper, we have confirmed 
that this is also true for ANM. 

Figure [H] shows the entropy per nucleon S/N of ANM derived from F/N: 
It can be seen that S/N increases with x monotonically. The approximate en- 
tropy per nucleon So/N is also shown: S/N and Sq/N are in good agreement 
with each other, too. 
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Figure 8: Entropies per nucleon, S/N, at T =10, 20 and 30 MeV for various values of x 
as functions of p. The approximate entropies per nucleon, Sq/N, are also shown. 



The pressure P of ANM is shown in Fig. It can be seen that P 
increases monotonically as x decreases. For a fixed x, P increases with T, 
and at T higher than a critical value Tq, P increases monotonically with 
p, i.e., at T < Tq, the mechanical instability occurs. In the case of SNM, 
To is equal to the critical temperature Tc of the liquid-gas phase transition 



32|, and Tq ~ 18 MeV, as reported in paper I. It is found that Tq decreases 



with X, and becomes MeV at x < 0.08. Here it should be noted that, in 
order to determine Tc of ANM, the chemical instability, in addition to the 



mechanical instability, has to be taken into account, as discussed in Ref. [32 
In our project to construct an SN-EOS, as mentioned before, we will treat 
this non-uniform mixed phase in the Thomas-Fermi approximation. 

Figure dO] shows the proton chemical potential fip and the neutron chem- 
ical potential /in at T = 20 MeV: In the case of SNM, /Xp = /in, and, as x 
decreases, /in decreases while /ip increases. It should be noted that /tp and 
/in shown in this figure are derived from F/N with the thermodynamic rela- 
tions, and are different from /toi introduced in Eq. f l25|) : The latter include 
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Figure 9: Pressures at T =10, 20 and 30 MeV for various values of x as functions of p. 



the effect of the single-particle potential energies that are not included in 
ei{k) given by Eq. (!26l) . 

The effective masses m* are shown in Fig. [TT]as functions of x. In the case 
of SNM, m* is identical to m*, and, as x decreases, m* decreases while m* 
increases. It is also seen that, as p increases, the T-dependence of m* becomes 
smaller. It should be noted that those m* are the variational parameters in 



the present variational method, but, as pointed out in Ref. [Ij], those can be 
regarded as weighted averages of the momentum-dependent effective masses 
of nucleons near the Fermi surfaces. In fact, the behavior of the obtained m* 
is fairly reasonable. 

Next, we discuss the validity of the frozen-correlation approximation 
(FCA) employed in this paper. This FCA is adopted in the original pa- 
per by SP, and the FHNC calculation in FCA is performed in Ref. |31l]- In 
principle, however, F/N should be minimized with respect to the correlation 
function fij as well as m*, as performed in the recent studies IJ, ll9|. There- 



fore, in the rest of this subsection, we perform the full minimization of F/N, 
i.e., for each set of m*, we minimize F/N with respect to fctsi^)^ ItA^) 
/gQ^(r) by solving the Euler-Lagrange (EL) equations derived from Eq. fl23|) . 
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Figure 10: The proton chemical potential /ip and the neutron chemical potential /in at 
T = 20 MeV for various values of x as functions of p. 
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Figure 11: Proton effective masses to* and neutron effective masses to* for T ~ 10, 20 
and 30 MeV at p = 0.16 and 0.32 fm~^ as functions of x. 
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Figure 12: Free energies per nucleon at T = 1 and 30 MeV for various values of x in the 
frozen-correlation approximation and those with the full minimization as functions of p. 



and then minimize F/N with respect to m*. In solving the EL equations, we 
impose two constraints as in the case at zero temperature, i.e., the extended 
Mayer's condition and the heahng-distance condition. The former is a con- 
dition similar to Eq. f lT^ . but Fp^^(r) and -FqFts('^) -^1^- dlSD ^^"^ 
replaced by F^^g{r]T) and F^p^^{r;T), respectively. In the latter condition, 
the healing distance is assumed to be the same as at zero temperature. 

In Fig. [121 the fully minimized values of F/N are compared with those 
of the FCA. It is found that the values of F/N in the FCA are in good 
agreement with those with the full minimization in a wide range of p, x and 
T. This agreement implies that the effects of the full minimization and the 
extended Mayer's condition cancel out each other. As a result, it is confirmed 
that the FCA is a good approximation in the present calculation, which is 



consistent with the result in Ref. 31 



3.3. Proton- fraction- dependence of the free energy 

In this subsection, we discuss the x-dependence of F/N. First, we define 
the following free energy difference: 

F(p,x,T) F(p,x = 1/2, T) , , 

x, T) ^ ^ - ^' ' ^ (29) 
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Figure 13: Free energy differences Fi/N given in Eq. (|29|) for various values of T at 
p = 0.16 and 0.32 fm"'^ as functions of 

and refer to Fsy^{p,T) = Fi{p,x = 0,T) as the symmetry free energy. At 
zero temperature, Fsym(p, T) reduces to the symmetry energy; i?sym = -^sym 
(po, T=0 MeV). 

Figure fT3l shows Fi{p,x,T) as functions of C = (1 — 2a;)^. Contrary to 
the case at zero temperature, it can be seen, in this figure, that Fi{p,x,T) 
are not proportional to (: The curves deviate from the straight hues more 
markedly as the temperature increases. 

For a more quantitative discussion, as in the case at zero temperature, 
we evaluate F{p, x,T)/N hj a linear interpolation between F/N of SNM and 
that of PNM with respect to (, i.e.. 



FUp,x,T) F(p,1/2,T) 



N 



N 



+ 



F(p,0,T) F(p,l/2,T) 



N 



(l-2x)2, (30) 



and then calculate the relative deviation AF{p, x, T) defined by 

^F(p...r)^ ^"'^---^^^_;^;''---^'/^ . (31) 

Figure [14] shows AF{p,x,T) as functions of (. As seen in this figure, the 
relative deviations at finite temperatures are much larger than those at zero 
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Figure 14: Relative deviations AF of the interpolated free energies from the original ones 
as functions of 



temperature. This implies that the quadratic interpolation between F/N of 
SNM and PNM with respect to x to obtain that of ANM, which is fairly 
good at zero temperature, is inappropriate at finite temperatures, as pointed 



out in Ref. 19 



In order to discuss the density dependence of the deviation mentioned 
above, we define the density-dependent symmetry free energy S'it(p, T) as 



^it(p,T) = 



1 F{p,x,T) 



(32) 

x=l/2 



It is obvious that Sit{p,T = OMeV) = S'i(p) at zero temperature. In Fig. 
[T5| we show Sit^PjT) as well as Fi{p,x,T) at T = 20MeV as a function 
of p. It can be seen that, in the wide range of p, Sn{p,T) is smaller than 
Fsyin(p, T) = Fi{p,x = 0,T) as in the case at zero temperature, but the 
difference is much larger, because the curvature of Fi{p, x, T) with respect to ( 
in Fig. [13] is much larger than that of Ei[p, x) shown in Fig. [31 Furthermore, 
since the thermal effect becomes smaller as p increases for a fixed T, the 
difference between Sit{p,T) and Fsym{p,T) decreases at high densities. 
Finally, in order to clarify the cause of the relatively-large deviation AF, 
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Figure 15: The density-dependent symmetry free energy Sit{p,T) defined by Eq. (j32l) at 
T = 20MeV as a function of the density p. Fi{p,x,T) a.t T — 20MeV are also shown for 
various values of x. 



we show, as functions of (, the following internal energy difference 

Eti(p, X, T) = , (33) 

and the entropy per nucleon S{p, x, T)/N, in Figs. [Tn]and[T7l respectively. As 
seen in Fig. Uni Et'i{p,x,T) is roughly proportional to (, while S{p,x,T)/N 
is not a linear function of ( as can be seen in Fig. [T71 Namely, the deviation 
of Fi from the straight line in Fig. [13] is caused by the x-dependence of S/N. 



4. Summary 

In this paper, we have constructed an equation of state (EOS) for uni- 
form asymmetric nuclear matter (ANM) at zero and finite temperatures using 
the variational method, which starts from the nuclear Hamiltonian with the 
AV18 and UIX potentials. At zero temperature, the obtained energy per nu- 
cleon E/N is reasonable as compared with that obtained by APR, and the 
mass-radius relation of NS calculated with the present EOS is consistent with 
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Figure 16: Internal energy differences E^i given in Eq. (j33p for various values of T at 
p = 0.16 and 0.32 fm"'^ as functions of 

observational data. We also confirmed the validity of the quadratic behavior 
oi E/N with respect to the proton fraction x at zero temperature. At finite 
temperature, the obtained free energy and other thermodynamic quantities 
such as internal energy, entropy, pressure and chemical potentials for ANM, 
which are important physical inputs in SN simulations, are reasonable. We 
also confirmed the thermodynamic self-consistency of the present calcula- 
tions, and the validity of the frozen-correlation approximation employed in 
this study. Finally, at finite temperatures, we have found a considerable de- 
viation of the free energy per nucleon F/N from the quadratic behavior with 
respect to x, and have shown that it is caused by the x-dependence of the 
entropy. 

In order to complete a new nuclear EOS table for SN simulations, the 
EOS of non-uniform nuclear matter is necessary. For this purpose, a Thomas- 
Fermi calculation for non-uniform matter, with use of the values of F/N for 
uniform matter constructed in this paper, is in progress. 
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Figure 17: Entropies per nucleon S/N for various values of T at p = 0.16 and 0.32 fm 
as functions of C. 



Appendix A. Explicit expressions for the functions employed in 
the energy formulas 

In this appendix, we present some of the explicit expressions used in this 
paper. At zero temperature, Fp^g(r), F^^p^^^r) and -F^Ftsl'^) ^'i^- &-^^M 
are given as 

isospinspin J 

= ^{60-(-l)*^%(ri2)/,(ri2)}, (A.l) 



K^ts{ri2) ^ E E Ul\L,2\'Pt'si2^Ydr,dr,---dr 

isospinspin J 



isospinspin 

2s + 1 f r 
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-(-1) 



1+8^12 



dr 



12 



dr 



12 



(A.2) 



FUs{n2) ^ E UULi2t PL2^Fdr,dr,---drN 

isospinspin J 



2s + 1 f r 



12 , 



70 
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12 



/i(ri2) 



+ /,(ri2) 



dHi{ri2) 
dr"^ 



12 



dr 



12 



^12 



(A.3) 



Here, ri2 = |''"i — ''"2I, ^ is the volume of the system and ^ represents 
the summation over the isospin and spin coordinates of all nucleons. The 
subscripts {i^j) in Eqs. ( 1A.ip -( tAl3|l represent (p,p), (p,n) and (n,n) for n = 
(+,0, — ), respectively. Furthermore, (Cp,^n) = — x) and 



h{r) = 3e 



jijkFtr) 



p,nj 



(A.4) 



At finite temperatures, the expression of E2T/N is similar to that of 
E2/N given by Eq. except that E^/N, F^Jr), F^^Jr) and F^^^Jr) 
defined above are replaced by En^/N, Fp^^{r;T), Fqpj^(r;T) and Fl^p^^{r;T) 
as follows: 



E- 



IT 



N 



2m TT^p 
2s + 1 



np{k)k dk + n^{k)k dk 
Jq 



4 



lf+%{r-T)l,{r;T) 



(A.5) 
(A.6) 



2s + 1 f 



4 



/ nj{k)k'^dk + / ni(k)k'^dk 



-(-1) 



fir 



(A.7) 
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+0 / ni{k)[k^ + —k^]dk 



+ 



3vr2p Jo 



3vr2p Jo 
dr 



nj{k)k'^dk 



1 

+ 2 



^ /,(r;T) dl,{r-T) ^ l^{r-T ) dk{r; T) 



dr 



dr 



Here, li{r;T) is defined as 



1 



TT^P Jo 



ni{k)jo{kr)Pdk {i = p, n). 



(A.8) 



(A.9) 



Appendix B. Validity of the interpolation by Lagaris and Pand- 
haripande for energies of asymmetric nuclear matter 

In this appendix, we examine the validity of the interpolation proposed 
by Lagaris and Pandharipande (LP) H for the calculation of E{p, x) /N in 
ANM at zero temperature; this interpolation is also adopted in paper 11. In 
this method, E{p,x)/N is approximated with use of E{p,x)/N for SNM and 
PNM as follows: 



^IntLp(p,x) _ Ei{p,x) ^ Ey{p, 1/2) ^ 



N 



N 



N 



Ev(p,0) Ey{p,l/2) 



N 



N 



'i-2xr. 

(B.l) 



Here, Ei(p, x)/N is the one-body kinetic energy given in Eq. dH]), and 

Ev(p,x) E{p,x) Ei{p,x) 



N 



N 



N 



(B.2) 



In Fig. m we show the relative deviation of EintLp{p,x)/N defined by an 
equation similar to Eq. (122|) . i.e.. 



AEi^p{p,x) 



Ej,,i^p{p,x)/N -E{p,x)/N 
Ei{p,x) 



(B.3) 
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It can be seen that \AEi^p{p,x)\ is much smaller than that of the usual 
interpolation Eq. fl2ip . because Ei{p,x)/N, which is the main cause of the 
relatively large values of AE, is removed from the interpolated part: The 
interpolation by LP is a good approximation for ANM at zero temperature. 
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